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The recent observations of the positron fraction in cosmic rays by PAMELA indi- 
cate that the fraction of positrons to the total electronic component in cosmic rays 
initially decreases in the energy region 1-10 GeV and increases thereafter. In this 
paper, we show that it is natural to expect such an increase of the positron fraction 
within the context of cosmic ray propagation models. It is shown that this ratio 
should reach an asymptotic value of ~0.6 at very high energies. The specific mea- 
surements by PAMELA help us to distinguish amongst various models for cosmic 
ray propagation, and in particular, they support the nested leaky box model. They 
also provide, in conjunction with the observations of the total electronic component 
by HESS, FERMI, ATIC, and other experiments, a way of estimating the spectrum 
of electrons directly accelerated by discrete sources of cosmic rays in the Galaxy. 



Direct observation of the cosmic ray electronic component (which includes both electrons 
and positrons with no charge discrimination) dates back to the early 1960s, and since that 
time, the energy range and the sensitivity of the observations have increased systematically. 
To date, we have at hand data from three new instruments, FERMI p], HESS [2]l3j, and 
ATIC [3], that have the requisite sensitivity to measure, with good statistical accuracy, the 
spectrum of the total electronic component (e~ +e + ) well into the TeV region. The reported 
spectrum in the region E >10 GeV is parameterized as 
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FIG. 1: Shown here is the total electron spectrum as measured by HESS, FERMI, ATIC, and other 
experiments. A smooth fit used for computational purposes is also shown. 

with T = 3.05 up to ~ 1 TeV, and the value of T increases to ~ 3.9 at higher energies. We 
reproduce in fig. [T] their observations and a compilation of the results of other measurements. 
The spectral slope below 10 GeV progressively flattens to a slope of ~1.7. Uncertainties in 
the fluxes are introduced due to solar modulation effects below an energy of a few GeV. We 
also show, in the same figure, a smooth fit to all the data that we adopt for some of the 
calculations. 

The electronic component in cosmic rays, because of its interactions with radiation fields 
such as starlight, the microwave background, and magnetic fields in the galaxy, has been 
particularly useful in understanding the origins and propagation of energetic particles in the 
Galaxy |5j. This and other early considerations of the effects on the spectral shape of the 
cosmic ray electrons were carried out in the context of a smooth distribution of cosmic ray 
sources in the galaxy, and the transport was described within the framework of the leaky 
box model [UJ. There were also attempts to calculate the flux of nuclear secondaries like 
Li, Be, and B, as well as positrons, in cosmic rays within the context of a nested leaky box 
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FIG. 2: The positron fraction as measured by PAMELA (+'s) is shown here. Also included is the 
calculation by Moskalenko and Strong [ID (solid line) and a smooth fit through the data (dashed 
line) . 

model which took into consideration the effects of storage of cosmic rays in a small bubble 
surrounding the compact sources of cosmic rays |7J [8 J . The key consideration for this model 
was the anisotropy of cosmic rays at high energies. If the residence time of cosmic rays in the 
Galaxy reduced with energy to accommodate the decreasing ratio of secondary to primary 
nuclei in cosmic rays, then the expected anisotropy would increase correspondingly at high 
energies. 

The recent observations of positrons in cosmic rays by PAMELA [9] has created much 
excitement because of the possible connection of these observations with annihilation of dark 
matter in the Galaxy. Their observations of the ratio R, of positrons to the total electronic 
component in cosmic rays is reproduced in fig. [2] The positron fraction at ~1.64 GeV was 
measured to be ~0.0673, which decreases to ~0.0483 at ~6.83 GeV and thereupon increases 
monotonically, reaching a value of 0.137 at a mean energy of 82.55 GeV. It is this monotonic 
increase that is being called anomalous, as it does not conform to the prediction of the 
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currently popular model of the cosmic ray propagation [10] [TT] . Accordingly, this paper 
begins with a review of the cosmic ray propagation models. The PAMELA measurements 
are to be viewed in the context of the measurements of the total electronic component and 
a compilation of the available data as shown in fig. [TJ 

The cosmic ray nucleonic component is dominated by protons with some neutrons coming 
in bound as He and other nuclei. The nucleon spectrum may be represented as 

f n (E) = K n E~P (2) 

where (3 ~ 2.6 — 2.7 (we adopt (3 = 2.65) in the energy region 1 GeV/nucleon to ~10 6 
GeV/nucleon beyond which the slope may increase to ~3. The p/n ratio effectively deter- 
mines the e + /e~ ratio generated by cosmic rays through interactions with matter in the 
sources and the interstellar medium through which they propagate before they leak out of 
the Galaxy. The theoretical calculation of e + /e~ generated through nuclear interaction of 
cosmic ray nuclei yield ~2 [12]. On the other hand, the observations of the ratio 
in cosmic rays gives ~1.3 [13]. We show our results for both of these values for the 

secondary positrons and electrons. 

In section II, we provide a brief overview of the models for cosmic ray propagation, and 
in section III, we discuss the PAMELA results, first in a model independent way, and then 
compare it with the expectations of the various models. Finally, section IV is devoted to a 
discussion of the main results of this paper and related matters. 

It should be noted that much work as has been done recently to explain the PAMELA 
data in terms of dark matter, pulsars, supernova remnants, and other astrophysical objects. 
Instead of listing the many papers which address these issues, we refer the reader to [H] 
which contains references to many recent works. 

II. BRIEF OVERVIEW OF MODELS OF COSMIC RAY PROPAGATION 

It is generally accepted that cosmic rays are accelerated in a large number of discrete 
sources distributed in the Galaxy, and the cosmic rays propagate from these sources moving 
along randomly oriented trajectories, akin to diffusion, until they leak away from the Galaxy. 
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During such a propagation, the cosmic rays might interact with the interstellar matter, the 
radiation, and magnetic fields. Any secondaries generated through such interactions, if 
charged, will be confined by the interstellar magnetic fields and will therefore follow the 
same kind of random paths as the primaries before escaping from the Galaxy. The various 
propagation models are characterized by the specific form chosen for the "vacuum path 
length distribution" [6J which describes the probability P(t) that the cosmic rays spend 
in any given region, such as a cocoon surrounding the sources, or in the general interstellar 
medium before escaping into the intergalactic space. The term "vacuum" emphasizes the fact 
that in specifying P(t), one considers hypothetical particles which do not suffer interactions 
or lose energy during propagation. The effects of these processes are to be added later on. 



A. The Leaky Box Model 

In its simplest original form [5J, one assumes that P(t) has a broad distribution with 
significant amplitude near t = 0, exemplified by a simple exponential function 

P(t) = e- t/T (3) 

where r is called the escape lifetime of the cosmic rays. In the original version, r was assumed 
to be sensibly independent of energy beyond ~l-2 GeV. Thereafter, since the discovery 
that the ratio of the fluxes of secondary cosmic ray nuclei to those of the primaries was a 
decreasing function of energy, r was considered to decrease with energy to accommodate 
the observations. We summarize in figs. [3] and [4] the available observations. The crucial 
issue here is how r behaves at energies beyond 10-20 GeV where the observations have low 
statistical significance or are non-existent at much higher energies. Most conventional models 
today [IP] [IT] assume that 

t a (E) ~ t E- a for E > 2GeV/n (Model A) (4) 



with A pa 0.4 — 0.5. Such an extrapolation to high energies is shown in a dot-dashed line in 
figs. |3]and[U 
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FIG. 3: The observed B/C secondary to primary ratio is plotted (points from a compilation in 
along with the power law extrapolation at high energies (dot-dashed line, Model A), a constant 
extrapolation (solid line, Model B), and a two-component fit (dotted lines, Model C). 

Alternatively, we may assume that t(E) becomes nearly constant at high energies 

t b {E) ~ t a (E) for E < WGeV 

tb{E) ~ tq ~ constant for E ^> IQGeV 

This is shown as a solid line in figs. [3] and |4| The latter form of r ~ tq at high energies 
predicts a lower level of anisotropy of cosmic rays compared to an ever decreasing ta(E). 
We will refer to the two models described briefly here as leaky box model A and leaky box 
model B respectively. 

B. The Nested Leaky Box Model 

An alternate way of accommodating the falling secondary to primary ratio is in the context 
of the nested leaky box model j8]. Here, one assumes that subsequent to acceleration, cosmic 
rays spend some time in a cocoon-like region surrounding the sources, interacting with the 
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FIG. 4: The observed (Sc+Ti+V)/Fe secondary to primary ratio is plotted (points from a compi- 
lation in along with the power law extrapolation at high energies (dot-dashed line, Model A), 
a constant extrapolation (solid line, Model B), and a two-component fit (dotted lines, Model C). 

matter there to generate some of the secondaries at lower energies. The residence time 
t s in the source region is energy dependant, decreasing with increasing energy. On the 
other hand, once these cosmic rays enter the general interstellar medium, their subsequent 
transport becomes independent of energy and the residence time becomes equal to tq. 

rJE) ~ r B (E) - Tq for 1 GeV < E < 10 GeV 1 
r G (E) ~ constant for 1 GeV < E < 10 6 GeV J 

The net effect of the interactions in these two regimes is to generate the correct ratio of the 
fluxes of secondary nuclei to those of their primaries. 

In this model, the anisotropics of cosmic ray fluxes remain constant and do not increase 
with energy. Moreover, the spectrum of cosmic ray primaries, say the nuclear component, 
is more easily understood. To see this, let S n (E) represent the injection rate of cosmic rays 
into the Galactic volume per unit volume and unit time at energy E per unit energy interval: 

S n (E) ~ E~ a . (7) 
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The spectrum f n (E) expected in the Galaxy is given by 



f n (E) ~ r G S n (E) ~ r G E 



—a 



(8) 



On the other hand, in the simple leaky box model A, 



f n (E) ~ 1 S n (E)r A ( J E) 
~ £-(«+ A ) 



for E>2 GeV/n 



(9) 



which will fit the observations for the choice a ~ 2.2 — 2.3 with the sources accelerating a 
flatter spectrum than the one that is observed. This flatter spectrum must continue up to 
very high energies, up to which the residence time continues to decrease as E~ A . If such a 
rapid decrease of residence time stops at any energy and becomes constant, then f n (E) will 
display E~ a behavior at higher energies, or until S n (E) itself changes its slope. 

In the leaky-box model B, the source function S n (E) should have an index a=2.2 at 
E < 10 GeV/n, which changes to a = 2.65 at E ^ 10 — 20 GeV, coincident with the change 
in behavior of tb{E), compensating its change and generating a smooth power law for f n (E). 
Thus, in the nested leaky box model, the observed spectral slopes simply correspond to that 
generated in the acceleration process in the cosmic ray sources. Also, with constant r G , the 
expected anisotropies do not increase with increasing energies, but remain sensibly constant. 

C. Spectrum of Secondary Electrons and Positrons in Cosmic Ray Propagation 



The generation of electrons and positrons in the interactions of the cosmic ray nuclear 
component occurs through the production of mesons, mainly pions, which decay to muons 
which in turn decay into electrons or positrons, transferring, on the average, a fraction of 
about 0.05 of the energy per nucleon of the primary. This is in contrast with the production 
of secondary nuclei, such as boron from the collision of carbon nuclei, where boron emerges 
with almost the same energy per nucleon as the primary carbon nucleus. This difference in 
their production characteristics leads to nearly identical source spectra S n - and S n + for the 
secondary electrons and positrons ~ E" 13 in all the three models: A, B, and C. 



Models 
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On the other hand, their equilibrium spectra f n + (E) and f n - (E) are markedly different 
in the three models. At energies where the energy losses due to synchrotron radiation and 
inverse Compton scattering on radiation fields are not important. The three models generate 
the spectra noted below: 

f n+ ~ S n+ (E)r A (E)^r E-^ +0 ^ (Model A) (fO) 

~ S n+ (E)T B (E) ~ r E-^ +0A) for E < 10 GeV (Model B) 

~ S n+ (E)T G ~ E' 13 for E > 10 GeV (Model B) 

or S n+ (E)r G ~ E^ 13 forE > 1 GeV (Model C) 



The spectra for the secondary electrons are similar to those given in eq. [10] except that 
because of the dominance of the protons in the cosmic ray beam, the production rate of 
positions is higher, with 

Sn - {E) = v . (11) 

This ratio r\ is theoretically estimated from the characteristics of high energy interactions to 
be ~0.5 [10J; on the other hand, the direct observation of the /i~/yU + ratio indicates a value 
of ~ 0.8 [T3]. In either case, rj is essentially independent of energy beyond a few GeV. 

Thus we see that, at high energies (E3>10 Gev), in the leaky-box model B and in the 
nested leaky box model C, the secondary positron and electron spectra are power laws with 
indices f3 equal to that of the spectrum of the nuclear component in cosmic rays. At very 
high energies, the energy losses due to synchrotron emission and inverse Compton scattering 
will steepen these spectra to f n + = E~^ +1 '. 

III. ANALYSIS OF THE POSITRON FRACTION OBSERVED BY PAMELA 

We find it useful to write the observed positron fraction R(E) in terms of the various 
components: 

R(E) = — fn+ . (12) 

Jn+ + In- + 9e- 

Here, f n + and f n - represent the positron and electron spectra generated as secondaries of 
the nuclear component of cosmic rays, and g e - is the spectrum of electrons resulting from 
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FIG. 5: Here we have subtracted f n +{E) and f n -{E) from the total spectrum of the electronic 
component and have shown g e (E), the spectrum of electrons generated by the cosmic ray sources. 
The primary electron component g e (E) is plotted for the HESS, FERMI, ATIC and other data 
using rj = 0.45. The positron spectrum from the nested leaky box model C is plotted as well (solid 
line) . 



direct acceleration in the sources. Note that there is no direct contribution to positrons from 
the source. It is convenient sometimes to work with the inverse of R(E) given by 

1 fn+ + fn~ + 9e- 



P{E) 



1 + 77 + 



9c 



R(E) f n + ' 

This allows one to find the spectrum of electrons generated by the sources g e - as 

g e -(E) = [P(E)-(l + r ] )]f n+ (E). 



(13) 



(14) 



This spectrum g e {E), generated in the Galaxy exclusively by the cosmic ray sources, is shown 
in fig. |5j 

Alternatively, we may just assume the functional form for f n +(E) given by the various 
propagation models and calculate the positron fraction by dividing this by ft(E), the total 
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FIG. 6: The theoretically calculated positron fraction in models A (similar to that of [TO] as shown 
in fig. [2]), B, and C are compared with the observations. All calculations are normalized at ~ 10 
GeV 



spectrum of electrons measured by FERMI, HESS, and other experiments, 

fn+(E) 



Rm(E) 



ft(E) 



(15) 



which is shown in fig. [6] along with the data from PAMELA. The normalization of the 
theoretical curves is such as to provide the best possible fit to the three models A, B, 
and C described earlier. (This normalization may indeed be explicitly calculated as it is 
proportional to ta(E), t b (E), and r G respectively for the three models and depends on the 
density of matter in the propagation region, the spectral flux of the nuclear component, the 
cross section for meson production, decay kinematics, etc.) In depicting the three model 
curves, we have included the effect of the energy losses at high energies. 

Comparison of the theoretical expectations of the positron fraction with the PAMELA 
data indicates that model A provides a rather poor fit to the observation, as already noted 
by several authors |9j. A careful calculation of the positron fraction under the general 
assumptions of model A was carried out over a decade ago by Moskalenko and Strong [TDJ. 
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Our estimates here are essentially the same as that derived by them. Even though both 
model B and model C predict nearly identical injection spectra, the equilibrium spectra at 
low energies differ drastically with each other. In model B, the injection spectrum has to 
be multiplied by r B (E) to get the equilibrium spectrum. On the other hand, in the nested 
leaky box model, we need to multiply only by r G [T] [HJ - Both these models predict identical 
equilibrium spectra at high energies, for E 3> 10 GeV. The good fit to observed the positron 
fraction shows that the residence time of cosmic rays is essentially independent of energy for 
E > 10 GeV, a constancy that is expected to continue up to ~ 10 5 GeV. 

In choosing between model B and model C, the latter is preferred from considerations 
of the spectra of primary nuclei as well. This is because for a simple power law input from 
the sources having a form S n (E) ~ E 13 , model B would be expected to yield a spectral form 
f n (E) ~ £-03+ A ) below ~ 10 GeV and E~ p at higher energies. On the other hand, the 
observed spectra of all the nucleonic components are simple power laws of slope ~ E~ 2 ' 65 
with no changes of slope in the tens of GeV region. Thus we conclude that the nested leaky 
box model provides good fit with the PAMELA data and is preferred also from consideration 
of other cosmic ray observations. 



IV. DISCUSSION AND RELATED MATTERS 



The main result that emerges from the present analysis is that the nested leaky box 
model provides a good fit to the positron fraction observed by PAMELA. The model is also 
consistent with other observations of cosmic rays. Until good measurements of the positron 
fraction was available, there was no easy means of choosing amongst the various models. The 
fact that the nuclear secondaries, such as Li, Be, and B, emerge from nuclear interactions 
with essentially the same energy per nucleon as their parents, C, N, and O, was the main 
cause for this uncertainty. However, the fact that the positrons carry, on the average, a 
fraction of only about 0.05 of the energy of their nuclear primaries breaks this degeneracy, 
allowing a choice to be made. Improvements of the measurements of the spectra of both 
the secondary nuclei and of positrons will help in fixing, more firmly, the parameters of the 
nested leaky box model. 
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The recent measurements of the electronic component of cosmic rays by the PAMELA, 
HESS, ATIC, and FERMI groups have opened up many interesting new avenues to discuss 
cosmic ray propagation. Some of these were discussed or hinted at in this paper. 
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